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Abstract - A class of exclusion processes in which particles perform history-dependent random 
walks is introduced, stimulated by dynamic phenomena in some biological and artificial systems. 
The particles locally interact with the underlying substrate by breaking and reforming lattice 
bonds. We determine the steady-state current on a ring, and find current-reversal as a function 
of particle density. This phenomenon is attributed to the non-local interaction between the walk- 
ers through their trails, which originates from strong correlations between the dynamics of the 
particles and the lattice. We rationalize our findings within an effective description in terms of 
quasi-particles which we call front barriers. Our analytical results are complemented by stochastic 
simulations. 



O ! Introduction. — Brownian motion can be rectified in 
— manifold ways [l}[4]- One design principle exploits inter- 
I actions between Brownian particles and a substrate lead- 
^> ing to local changes in the substrate properties which in 
GO turn affect the particle's motion. For example, the pro- 
cessive uni-directional motion of collagenase on collagen 
^} fibrils [5j[6] has been argued to be due to a burnt-bridge 
.mechanism [7]: collagenase, as it diffuses along collagen 
fibrils, cleaves its track at some recognition sites so that 
^"Tit always ends up behind the cleavage site. This then 
l acts as a burnt bridge biasing collagenase motion. Simi- 
r* larly, motivated by the translocation of Holliday junctions, 
• £h asymmetric nucleation of hydrolysis waves has been pro- 
posed as a track-mediated mechanism driving directed mo- 
5-ntion [8j|9]. Recently, artificial systems employing bipedal 
DNA motors have been engineered using a track-mediated 
rectification principle [10 . Moreover, synthetic molecu- 
lar systems, termed molecular spiders 
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have been 

constructed utilizing a similar mechanism. They consist 
of an inert body and catalytic legs. The legs are DNA 
enzymes which, upon binding to a complementary DNA 
substrate, cleave it into two shorter products that then 
have a lower affinity for the legs. The spider's motion 
thus changes the biochemical properties of the molecular 
track, and its trail consists of sites with enhanced dis- 



sociation rates. When interacting with precisely defined 
environments such molecular spiders even perform some 
elementary robotic behavior 
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Both of these examples are non-Markovian stochastic 
processes: through cleaving the track or changing kinetic 
parameters of the substrate the system acquires a mem- 
ory of the path traced out by the particle. Such pro- 
cesses belong to a general class of random walks on lattices 
with weighted bonds or sites. Well-studied models include 
the self-avoiding random walk [15], the reinforced ran- 
dom walk 
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and the excited random walk 
random walks with memory show unusual behavior, like 
anomalous diffusion or spatial confinement. 

We are interested in collective phenomena emerging 
from the interaction of many non-Markovian random 
walkers. Such systems are genuinely distinct from their 
Markovian counterparts like, e.g., the symmetric or asym- 
metric exclusion processes (ASEP) [18| or Brownian ratch- 
ets [19] . Whereas in the latter case interactions are local in 
space and time, they are non-local for non-Markovian ran- 
dom walkers because they mutually influence each other 
through their trails, i.e., changes induced in the track while 
moving over it. Since this implies strong correlations be- 
tween the dynamics of the particles and the lattice, we 
expect novel behavior as compared to Markovian systems 
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where inhomogeneities are either quenched, i.e., static on 
the time scale of particle dynamics, or annealed, i.e., fast 
and decoupled from the particle dynamics. We focus on 
the asymptotic dynamics at long times, where a stationary 
state is reached in both the particle and lattice dynamics. 
The main quantity of interest is the current as a function 
of the model parameters such as the kinetic rates and the 
particle density. We find a particularly interesting phe- 
nomenon: a current reversal with respect to the particle 
concentration. 

The model. — In this work, we generalize the re- 
versible "Burnt-Bridge model" (BBM) [7| [20}|23] to a non- 
Markovian many body system where particles interact 
through their trails and on-site exclusion. In analogy to 
the asymmetric simple exclusion process (ASEP) we term 
the model burnt-bridge exclusion process (BBEP). In this 
model a fixed number N of random walkers moves along 
a periodic one-dimensional lattice with sites i — 1, . . . , L 
connected by bonds which can be in either of two states: 
intact or broken; see fig. [I] Each site is either occupied by 





Fig. 1: Illustration of the burnt-bridge exclusion process. The 
one-dimensional lattice consists of bonds, which can be either 
in an intact (dashed green) or in a burnt state (zig-zag red). 
Only upon crossing an intact bridge to the right it gets burnt. 
Passing a burnt bridge to the right or left repairs a bridge at 
rate r±. Particles are subject to on-site exclusion. 

a single particle or is empty. Bonds are subject to local 
interaction with crossing particles which may "destroy" a 
formerly intact bond, leaving it in a broken state ("burnt 
bridge") or "heal" a burnt bridge restoring it into an in- 
tact state again. If a bond is intact a particle may cross 
it with rate 1 in either direction (fixing the time scale). 
We assume a totally asymmetric burning mechanising] A 
particle will always burn an intact bond when crossing it 
in the positive direction (to the right), while it will never 
damage an intact bond, when crossing it in the negative di- 
rection (to the left). Once a bridge has been burnt, it poses 
an obstacle for any approaching particle. However, these 
are only transient obstacles for the particle traffic: We as- 
sume that broken bonds are "repaired" at reduced rates 
r± < 1, that depend on the crossing direction as explained 
in fig. [T] For the parameter regime r_ ^> r + the asym- 
metric repair mechanism and the burnt-bridge mechanism 
steer the particles in opposite directions. In this case, the 




Fig. 2: Typical trajectory of a random walker (blue) with 
repair rates r_ = 0.65, r+ = 0.15, and lattice size L — 10. 
Horizontal lines indicate the state of the bonds: intact (dashed 
green) or burnt (solid red). The "front barrier" is marked by a 
thick red line. Pushing events occur at times to,...,ts, and two 
ling events happen right after t = ta. 



non-Markovian character of the random walk may show 
its most explicit manifestation as a current reversal, both 
with respect to the lattice length and the number of par- 
ticles. 

The front barrier picture. — Already the dynam- 
ics of a single walker pose a highly difficult mathematical 
problem which can be solved exactly only in some special 
Here we show that for the parameter regime, 
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1 This assumption is made for simplicity. All main conclusions of 
our analysis remain valid for more complicated versions of the model 
[241. 



cases 

r_ 3> r + , where the asymmetric repair mechanism and the 
burnt-bridge mechanism are antagonists, there is an effec- 
tive description in terms of a kind of quasi-particle which 
we call the "front barrier". To develop this picture we 
start our analysis by exploring the characteristic features 
of a random walker's trajectory; see fig. [2j We start the 
particle's journey on a lattice with all bonds broken such 
that it finds itself trapped between two partially reflecting 
barriers. Since r_ 3> r + broken bonds are more "transpar- 
ent" in the backward direction, the walker most probably 
bounces off the broken bond in the forward direction, and 
repairs the bond in the backward direction, c.f. fig. [2] 
We call the strongly reflecting bridge in the forward direc- 
tion "front barrier" (FB). In every step backward bridges 
are repaired, while they are burnt for each forward step. 
Thus, the particle finds itself not only initially, but for 
some extended time period in front of a burnt bridge and 
behind the FB; here during [t , ti] . Since r_ < 1 the ran- 
dom walker is biased to move towards the FB but bounces 
off several times before it finally happens to repair it at 
time t\. Now, it finds itself in a similar situation as at 
time io; he., behind a strongly reflecting barrier. With 
an intact bridge behind, it most probably again bounces 
off the reflective wall. When doing so, the particle does 
not burn the intact bridge, as this is a backward crossing 
event. Up to this minor detail the particle dynamics dur- 
ing the time interval [£1,^2] follows the same principles as 
for [toj^i]- The FB has now been pushed forward by the 
walker repairing it at t\. 

This pattern of movement could basically repeat itself 
indefinitely; see time intervals [toi*i]) [ii>i2]>---> [^4^5] in 
fig. [2] However, since the lattice is finite, the random 
walker may at some point in time also travel the whole 
length of the ring in the backward direction, and eventu- 
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Fig. 3: (a) Front barrier picture as explained in the text, (b) 
The velocity v of a single random walker as a function of lattice 
length L. Simulation results (symbols) are given for different 
sets of repair rates. Solid lines are obtained from the front 
barrier picture, eq. ph. 



ally arrive on the weakly reflecting side of the FB. From 
there, it can easily pierce (and thereby remove) the FB. 
At this time, t — t$, there are no more burnt bridges left 
on the lattice, and the random walk becomes symmetric, 
but only until crossing of a bridge in the forward direc- 
tion induces this bridge to burn. In our sample trajectory 
the random walker makes two backward steps before it 
turns forward again. It leaves behind a trace of burnt 
bridges and eventually arrives at a newly formed, strongly 
reflecting FB, which has been pulled back by two steps. In 
summary, the dynamics of the FB is determined by the an- 
tagonism between pushing and pulling events. This leads 
to a current reversal as a function of the lattice length 
L; cf. fig. [3j For long rings, L> 1, the particle, due to 
the weak bias in the positive direction, is localized in the 
immediate vicinity of the front barrier. Then, the dynam- 
ics are dominated by pushing events, and the velocity is 
positive. In contrast, shortening the lattice increases the 
relative frequency of pulling events. The probability of 
such events heavily depends on the random walker's abil- 
ity to diffuse over the whole length of the lattice. There- 
fore, for small enough lattices, pulling events may at some 
point outnumber pushing events, resulting in a negative 
velocity. 

These heuristic observations can also be formulated 
quantitatively. For the parameter regime r + <C r_ < 1 
the state space most relevant for the dynamics is highly 
reduced. Instead of all the 2 L possible bridge configura- 
tions, only "block configurations" are statistically signifi- 
cant. These consist of two continuous regions on the lat- 
tice where the bridges are either burnt or intact, a "burnt 
block" and an "intact block", separated by a front bar- 
rier. The main reason for this reduction is a separation of 
time scales: since repair of bridges in the forward direc- 
tion are rare events (r+ <C r_ < 1) a random walker finds 
itself most of the time behind a strongly reflecting bridge, 
the FB^] In this FB picture the random walker's effective 



2 For large r+, the approximation fails since it becomes more and 
more likely that a random walker pushes the FB forward by not only 
one but many steps. Then, after reflection at the FB, it is able to 
constitute a new block of burnt bridges inside the previously intact 
block, and the picture with a single block configuration breaks down. 



states are a = 1, . . . ,L, where it is located at a distance 
a behind the FB, and a = L + 1, where it has crossed the 
FB from the right and finds itself on a fully intact lattice; 
see fig. [2] and fig. [3^,. In the latter state the FB is pulled 
backward at rate 1, while in state a = 1 the FB is pushed 
forward at rate r+. Hence the velocity of the FB and with 
it the average particle velocity is given by 



v = r + ■ mi - 1 • m L+ i , 



(1) 



where m a denotes the probability to find the walker in 
state a. Since the boundary states a = 1, L + 1 are to- 
tally reflecting, the probability current between two states 
is strictly zero: j a = 1 ■ m a+ i — r_ ■ m a — 0. This im- 
plies that m a decays exponentially, m a = {r-) a 1 mi. 
With the normalization of the probability distribution, 
J2a=l m a = 1, this gives mi = (1 — r_)/(l — r^ +1 ) such 
that the average particle velocity eq. 0, reads 



v(L) 



l-r 



L+l 



(2) 



These results can be generalized to lattices containing 
both strong and weak links, where strong links are un- 
affected by the random walker and weak links correspond 
to the bridges considered here [24] . 

Equation |2]) quantifies the relative frequency of FB 
pushing and pulling events, and is compared with results 
obtained from stochastic simulations of the BBM in fig.[3]3. 
As expected there is nice agreement for r+ <C r_ < 1; the 
analytical result is even exact for r + = 0. This allows us 
to use the FB picture, eq. ([2]), to determine the precise 
conditions for velocity reversal: 



< r + < r_ < 1 . 



(3) 



In this regime the velocity becomes zero at a length 
Lq = In r + /In r_ , which may be interpreted as a balance of 
forces exerted on the particle by a broken bond to its front 
and rear, respectively. Both of these results agree very 
well with our numerical results, summarized as a phase 
diagram in fig. [4^,. 

The burnt-bridge exclusion process. We now 

turn to investigate collective properties emerging when 
there are many random walkers present at a finite den- 
sity p. These interact in two ways: local on-site exclusion 
and non-local interaction through the history of the ran- 
dom walkers, i.e., the sequence of intact and burnt bridges 
each walker leaves behind in its trail. The effect of on-site 
interaction alone on the collective properties is well under- 
stood. In the ASEP with periodic boundary conditions, 
the density p is homogeneous and there are no correlations 
in the steady state 18 . Hence the current follows from 



a mean-field argument as the product of the difference in 
forward and backward hopping rates, A±, and the proba- 
bility to find a given site occupied and its neighboring site 
empty: j = (A + — A_)p(l — p). The sign of the current 
only depends on the difference between the forward and 
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Fig. 4: Phase diagrams for (a) a single random walker on a pe- 
riodic lattice with varying length L and (b) varying number of 
random walkers and fixed L = 160. Dashed lines represent the 
phase boundary calculated from the FB picture, eq. d3b. Solid 
line is a guide to the eye indicating the actual phase bound- 
ary. Employing a Gillespie algorithm, the parameter regime 
has been scanned with a resolution of 8r± — 0.05. Red crosses 
(+) indicate a phase with positive velocity, blue crosses (x) a 
phase with negative velocity, and green stars (*) a phase with 
current reversal. 



backward rate, but not on the particle density. The latter 
regulates the magnitude of the current. 

In the BBEP with both on-site exclusion and trail- 
mediated interaction between the random walkers the be- 
havior is much richer. In fig. [5] simulation results for the 
average velocity v(p) :— j/p as a function of the particle 
density p are shown for both the ASEP and the BBEP. 
For r + 3> r_ , we observe a roughly linear velocity-density 
relation, v cx 1 — p, closely resembling the relation for 
the ASEP but with a reduced amplitude. In contrast, for 
r + <C r_ , we find a nonlinear dependence of the velocity 
on the density p strongly deviating from the correspond- 
ing ASEP result which would give a negative velocity. In 
particular, for some parameter ranges, there is current re- 
versal at two critical values of the particle density. 

We attribute the nonlinear v (p)-relation for r + -C r_ to 
the history-dependence of the random walks. Every ran- 
dom walker changes the weights of the lattice by either 
repairing or burning bridges as it traverses the lattice. 
This affects the dynamics of neighboring particles walk- 
ing on its trail. What really matters is, however, only the 
fact that the weights, i.e., state of the bridges, have been 
changed but not when these changes had occured. This 
suggests that a walker reacting to the trail of its nearest 
neighbors might behave quite similar to a single random 
walker placed on a periodic lattice, always encountering 
its own trace. In other words, the dynamics of random 
walkers traveling at a mean distance d should strongly re- 
semble those of a single particle on a ring of length d, with 
identical repair rates r±. 

To illustrate this argument, we compare trajectories of 
random walkers for a few and many particles in fig. [6j 
With each particle we can, at any point in time, asso- 
ciate an unambiguously defined FB. With two random 
walkers on the lattice, direct encounters are observed only 
rarely. Most of the time, both are continuously pushing 
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Fig. 5: Veloctty-density relation, (a) Average particle velocity 
for the BBEP and ASEP as a function of particle density p 
for a set of kinetic parameters indicated in the graph; L = 
40. Simulations were performed using the Gillespie algorithm, 
(b) Comparison of simulations with the effective single-particle 
approximation (solid lines) for r_ = 0.8 and r+ = 0.2. 



their own FB forward without ever taking notice of the 
other walker's presence. Only in one rare event, the upper 
walker reaches the FB of its companion and pulls it back 
by one lattice unit. This happens only occasionally, and 
for large times, the average velocity is positive (in the for- 
ward direction). In contrast, if there are many particles 
on the lattice pulling events become more frequent sim- 
ply because there is a shorter mean distance between the 
walkers. As a consequence, for the parameter values cho- 
sen in fig. [6j the average particle velocity is reversed since 
pulling events outnumber pushing events of FBs. 

Upon replacing the lattice length L by the mean particle 
spacing d — 1/ pin eq. |2| one obtains an effective single- 
particle approximation for the many-body system. As il- 
lustrated in fig.[5j3, it gives a reasonable description of the 
velocity-density relation at small densities but becomes 
quantitatively wrong with increasing density. This is to 
be expected since correlations due to on-site exclusion are 
not accounted for in the FB picture. The results obtained 
from the FB picture for low densities can be generalized to 
very high densities close to p = 1 upon employing particle- 
hole symmetry, which implies pv(p) = (1 — p)v(l — p). 

The single-particle front barrier picture qualitatively ex- 
plains the observed current reversal in the many-particle 
BBEP. It identifies competition between pushing and 
pulling events of the FBs as the main mechanism. This 
trail-mediated mechanism is distinct from other current 
reversal mechanisms observed in Brownian ratchet models 
19 . There an interplay of the local ratcheting mechanism 
and hard-core repulsion between the particles leads to a 
reversal of the net current with increasing particle density. 
While in the latter case, excluded volume effects are driv- 
ing current reversal, they actually weaken it in the BBEP. 
As can be inferred from fig. [5j) the actual reduction in the 
average particle current with density is less than the FB 
picture would predict. Crowding effects asymmetrically 
affect the frequency of pulling and pushing events, to the 
disadvantage of pulling events. 

To map out the phase diagram we performed exten- 
sive stochastic simulations for the BBEP and measured 
velocity-density relations for a broad range of repair rates; 
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Fig. 6: Typical particle trajectories of the BBEP with repair 
rates r_ = 0.85, r+ = 0.05, and lattice size L = 60 for two 
(top) and six (bottom) particles. The trajectories are alter- 
nately drawn in blue (dark grey) and turquoise (light gray). 
For better visualization of the trajectories the bridges' states 
are not shown. Front barriers are only shown in the top but 
not in the bottom figure. 



see fig. [4)3. Indeed, there is a broad parameter regime, 
where current reversal is observed; crowding effects, how- 
ever, narrow this regime as compared to the results ob- 
tained from the effective single-particle picture which con- 
siders trail-mediated effects only. 

Summary. — The burnt-bridge exclusion process has 
features which are genuinely distinct from the Markovian 
analog asymmetric exclusion process. The most promi- 
nent is current reversal as a function of density. This phe- 
nomenon is absent for the asymmetric exclusion process 
and is attributed to the non-local interaction between the 
walkers through their trails. Those interactions are caused 
by a strong coupling between the dynamics of the parti- 
cles and the lattice which lead to memory effects acting 
opposite to the broken bonds' asymmetric transmissibility. 
The observed phenomenon is robust to model variations 



24 . We suppose that exclusion processes for the broader 
class of random walks on lattices with weighted bonds or 
sites shows an even richer dynamics with a plethora of new 
phenomena, e.g. by considering non-Markovian analogs of 
generalized exclusion processes 25-27 . It will be also in- 
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